Nucleon electromagnetic form factors in twisted mass lattice QCD 
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We present results on the nucleon electromagnetic form factors within lattice QCD using two 
flavors of degenerate twisted mass fermions. Volume effects are examined using simulations at two 
volumes of spatial length L = 2.1 fm and L = 2.8 fm. Cut-off effects are investigated using three 
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extrapolated to the physical pion mass allowing for a comparison with experiment. 
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I. INTRODUCTION 



Understanding the structure of the nucleon using the 
underlying theory of the strong interactions is a central 
problem of hadronic physics. The nucleon (N) electro- 
magnetic form factors provide an indispensable probe to 
the structure of the nucleon. Experiments to measure 
the electromagnetic nucleon form factors have been car- 
ried out since the 50's. A new generation of experiments 
using polarized beams revealed unexpected results [l|, . 
The form factors obtained in these polarization exper- 
iments differ from those extracted in previous experi- 
ments based on the Rosenbluth cross-section separation 
method. The new generation of experiments has shown 
that the ratio of the proton electric to magnetic form fac- 
tor G^/G^j decreases almost linearly with increasing mo- 
mentum transfer squared instead of being approximately 
constant. Two-photon exchange effects, previously ne- 
glected, were shown to be the source of the discrepancy. 
For a recent review we refer the reader to Ref. 0, 0] . Pre- 
cision experiments are currently under way at major fa- 
cilities in order to measure the nucleon form factors even 
more accurately and at higher values of the momentum 
transfer 

In this work we present results on electromagnetic 
form factors obtained using two degenerate light quarks 
{Np—2) in the twisted mass formulation. Twisted mass 
fermions (TMF) d Q provide an attractive formulation 
of lattice QCD that allows automatic 0{a) improvement, 
infrared regularization of small eigenvalues and fast dy- 
namical simulations For the calculation of the nu- 
cleon form factors, which is the aim of this work, the au- 
tomatic 0{a) improvement is particularly relevant since 
it is achieved by tuning only one parameter in the action. 



requiring no further improvements on the operator level. 

The action for two degenerate flavors of quarks in 
twisted mass QCD is in lattice units given by 

S ^ Sg + ^ X{x) [Dw+mci-it+ij5T^ fj] xi^) , (1) 



where Dw is the Wilson Dirac operator. For the gluon 
sector we use the tree-level Symanzik improved gauge 
action, Sg Q. The quark fields x ^"^^ in the so-called 
"twisted basis" obtained from the "physical basis" at 
maximal twist by a simple transformation: 
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and ^^x^i^ + ir^ls] . (2) 



We note that, in the continuum, this action is equivalent 
to the standard QCD action. A crucial advantage is the 
fact that by tuning a single parameter, namely the bare 
untwisted quark mass to its critical value Wcrit, a wide 
class of physical observables are automatically 0{a) im- 
proved. A disadvantage is the explicit flavor symmetry 
breaking. In a recent paper we have checked that this 
breaking is small for the baryon observables under con- 
sideration in this work and for the lattice spacings that we 
use Simulations including a dynamical strange 

quark are also available within the twisted mass formula- 
tion. Comparison of the nucleon mass obtained with two 
dynamical flavors and the nucleon mass including a dy- 
namical strange quark has shown negligible dependence 
on the dynamical strange quark [l^. We therefore ex- 
pect the results on the nucleon form factors to show little 
sensitivity to a dynamical strange quark as well. 

In this work we consider simulations at three values 
of the coupling constant spanning lattice spacings from 
about 0.05 fm to 0.09 fm. This enables us to examine 
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the continuum limit of the electromagnetic form factors. 
We find that cut-off effects are small for this range of 
lattice spacings. We also examine finite size effects by 
comparing results on two lattices of spatial length L = 
2.1 fm and L = 2.8 fm [iMl. 
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II. LATTICE EVALUATION 



Correlation functions 



FIG. 1: Connected nucleon three-point function. 



To extract the nucleon form factors we need to eval- 
uate the nucleon matrix element {N{p' , s')\jfj_\N{p, s)) , 
where \N(p',s')), \N{p,s)) are nucleon states with final 
momentum p' and spin s', and initial momentum p and 
spin s. 

The nucleon electromagnetic matrix element for real 
or virtual photons can be written in the form 



N{p',s') \r\N{p,s)) 



I'M 



ENip') ENip) 



1/2 



u{p',s')0''u{p,s) 



(3) 



where — (p'—p)'^, rriM is the nucleon's mass and i?jv(p) 
its energy. 

The operator can be decomposed in terms of the 
Dirac and Pauli form factors as 



2tojv 



(4) 



where i^i(O) = 1 for the proton and zero for the neutron 
since we have a conserved current. ^2(0) measures the 
anomalous magnetic moment. They are connected to the 
electric, Ge, and magnetic, Gm, Sachs form factors by 
the relations 



Guiq") 



(5) 



An interpolating field for the proton in the physical 
basis is given by 



J{x) 



[u'^^(x)C75d''(x)] u%x) 



(6) 



and can be written in the twisted basis at maximal twist 
as 
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abc 



u"' {x)Cj5d''{x) rix). (7) 



The third component of the isovector current is invariant 
under rotation from the physical to the twisted basis. 

In order to increase the overlap with the proton state 
and decrease overlap with excited states we use Gaussian 
smeared quark fields [l^ [l^ for the construction of the 



interpolating fields: 

C.car(i, x) = ^ F^'i^, y; U{t)) q\t, y) , (8) 
y 

F = + 
3 

4=1 

In addition, we apply APE-smearing to the gauge fields 
Ufj, entering the hopping matrix H. The smearing pa- 
rameters are the same as those used for our calculation 
of baryon masses with a and n optimized for the nucleon 
ground state [l^. The values are: a = 4.0 and n = 50, 
70 and 90 for ^ = 3.9, 4.05 and 4.2 respectively. 

In order to calculate the nucleon matrix element of 
Eq. ([3]) we calculate the two-point and three-point func- 
tions defined by 

G{ci,tf)=Y^ e-^-r^T^^- {J^{tf,^f)J^{t,,^,)) (9) 

where Fq and F^ are the projection matrices: 

ro = ^(l+7o), Ffc=iro757fc- (11) 

The kinematical setup that we used is illustrated in 
Fig. [T] The creation (source) operator at time ti=0 has 
fixed spatial position Xi=0. The annihilation (sink) op- 
erator at a later time tf carries momentum p'=0. The 
current couples to a quark at an intermediate time t and 
carries the momentum q. Translation invariance enforces 
q = — p for our kinematics. The form factors are calcu- 
lated as a function of — —q^ > 0, which is the Eu- 
clidean momentum transfer squared. Provided the Eu- 
clidean times, t — ti and tf —tt are large enough to filter 
the nucleon ground state, the time dependence of the Eu- 
clidean time evolution and the overlap factors cancel in 
the ratio 



i?'^(F,q,t) 



G^(r,q,t) G{p,tf~t)G{0,t)G{0,tf) 
GiO,tf) Y G(0,t/-i)G(p,t)G(p,t/) ' 

(12) 
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yielding a time-independent value 

lim lim i?^(r,q,i) =n^(r,q). (13) 

We refer to the range of i-values where this asymptotic 
behavior is observed within our statistical precision as 
the plateau range. We use the lattice conserved electro- 
magnetic current [33}. j''(x), symmetrized on site x by 
taking 
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(14) 



with an additional set of sequential inversions. Since in 
this work, we consider the temporal and spatial F's we 
need a total of four sets of sequential inversions. 

The nucleon matrix element also contains isoscalar vec- 
tor current contributions. This means that disconnected 
loop diagrams also contribute. These are generally dif- 
ficult to evaluate accurately, since the all-to-all quark 
propagator is required and the signal to noise ratio is 
extremely low. In order to avoid disconnected diagrams, 
we calculate the isovector form factors. Assuming SU (2) 
isospin symmetry, which holds to 0{a?) in the twisted 
mass formulation, it follows that 



We can extract the two Sachs form factors from the ratio 
of Eq. by choosing appropriate combinations of the 
direction /i of the electromagnetic current and projection 
matrices F. 

Inclusion of a complete set of hadronic states in the 
two- and three-point functions leads to the following ex- 
pressions, written in Euclidean time: 
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n'^='(rfe,q) = C- e,,fe Gm{Q 



(15) 



1 - 2 _ 1 - 

-d'y''d)\p) - { 7i|(-M7^M - -d7''d)|n) 



(19) 



One can therefore calculate directly the three-point func- 
tion related to the right hand side of the above relation 
which provides the isovector nucleon form factors 



GliQ^), 
GUQ')- 



(20) 



n'^=^(ro,q) = c 



2m N 



Gb(Q2) 



n^=°(ro,q) = c 



En + mN 
2mN 



GsiQ^) 



(16) 



(17) 



where C 



2mi 



is a kinematical factor con- 



-Eiv(-Eiv+miv) 

nected to the normalization of the lattice states and the 
two-point functions entering in the ratio of Eq. (|12p [2l| . 
The first observation regarding these expressions is that 
the polarized matrix element given in Eq. (jisp . from 
which the magnetic form factor is determined, does not 
contribute for all momenta q. New inversions are neces- 
sary every time a different choice of the projection matrix 
Ta is made and therefore to get the other components we 
would need two additional inversions. Alternatively, one 
can construct a suitable linear combination for the nu- 
cleon sink that leads to [211 



nSJta) = i:n»-(r.,,) = ^{ 

fe=l 



(93 - 92)<5i,j 



+ {qi - qz)52., + (92 - <Zi)53,«|GM(g') (18) 

which is optimal in the sense that it provides the maxi- 
mal set of lattice measurements from which Gm can be 
extracted, requiring one set of sequential inversions. One 
can choose the sink of Eq. (ITS)) or do three inversions one 
for each spatial F^. Which choice is more cost effective 
needs to be determined by comparing the statistical error 
at fixed cost. For the evaluation of the electromagnetic 
form factors the two options are almost equivalent. No 
such improvement is necessary for the unpolarized ma- 
trix elements given in Eqs. ([TO]) and PT)). which yield Ge 



The isovector electric form factor, G^^", can be obtained 
directly from the connected diagram shown in Fig. [T] 
To extract this quantity we consider either the spatial 
components of the electromagnetic current as given in 
Eq. (|16p or the temporal component given in Eq. (|17|) . 
The isovector magnetic form factor, " is extracted 
using Eq. ()15|1 for all three spatial components. 

Besides using an optimal nucleon source, the other im- 
portant ingredient in the extraction of the form factors is 
to take all the lattice momentum vectors that contribute 
to a given into account in our analysis. 

If a form factor G{Q^) can be extracted according to 
Eq. (fT5|) - P?|) from a total of M directions fi and lattice 
momenta q and we denote the plateau values by Ilfc , their 
statistical errors by Wk and the corresponding coefficient 
by Gfc, the form factor is calculated by minimizing 



M 

El 

fe=i 



CfcG(Q^) - n, 

Wk 



(21) 



This is a least-squares fit to a constant and the result is 
the weighted average of the individual measurements. 

Collecting contributions from all q directions im- 
proves the statistical precision and is moreover necessary 
to guarantee automatic 0(a)-improvement with twisted 
mass fermions. Phenomenologically interesting quanti- 
ties like the r.m.s. radii and magnetic moments can thus 
be obtained with increased precision. 

The connected diagram Fig. [T] is calculated by per- 
forming sequential inversions through the sink yielding 
the form factors at all possible momentum transfers and 
current orientations /i. Since we use a sequential inver- 
sion through the sink we need to fix the sink-source sep- 
aration. Statistical errors increase rapidly as we increase 
the sink-source separation. Therefore we need to choose 
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the smallest possible that still ensures that the nucleon 
ground state dominates when measurements are made 
at values of t in the plateau region. In order to check 
that a sink-source time separation of '-^ 1 fm is sufficient 
for the isolation of the nucleon ground state we compare 
the results at /3=3.9 obtained with {tf — ti)/ a— 12 i.e. 
{tf — ti) ^ 1 fm with those obtained when we increase 
to {tf — ti)/a=14: [l^H^. It was demonstrated that the 
plateau values for these two time separation are com- 
patible yielding the same results. This means that the 
shorter sink-source separation is sufficient and the ground 
state of the nucleon dominates in the plateau region. We 
therefore use in all of our analysis tf — ti ^ 1 fm. 

B. Simulation details 

The input parameters of the calculation, namely /?, 
L/a and a/i are summarized in Table HI The lattice spac- 

I 



13 = 3.9, a = 0.089(1)(5) fm, ro/a = 5.22(2) 


24^ X 48, L = 2.1 fm 
32^ X 64, L = 2.8 fm 


0.0040 0.0064 0.0085 0.010 
Stat. 944 210 365 477 

(GeV) 0.3032(16) 0.3770(9) 0.4319(12) 0.4675(12) 
m^L 3.27 4.06 4.66 5.04 
afi 0.003 0.004 
Stat. 667 351 

(GeV) 0.2600(9) 0.2978(6) 
m^L 3.74 4.28 


13 = 4.05, a = 0.070(1)(4) fm, ro/a = 6.61(3) 


32^ X 64, L = 2.13 fm 


a/i 0.0030 0.0060 0.0080 
Stat. 447 325 419 

(GeV) 0.2925(18) 0.4035(18) 0.4653(15) 
m^L 3.32 4.58 5.28 


P = 4.2, a = 0.056(1)(4) fm ro/a = 8.31 


32^ X 64, L = 2.39 fm 
48^ X 96, L = 2.39 fm 


an 0.0065 
Stat. 357 

(GeV) 0.4698(18) 
m^L 4.24 
a^i 0.002 
Stat. 245 

(GeV) 0.2622(11) 
m-nL 3.55 



TABLE I: Input parameters {/3,L,an) of om lattice calculation and corresponding lattice spacing (a) and pion mass (m-^). 
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C. Determination of the lattice spacing physical units. The nucleon mass has been computed on 

the same ensembles that are now used here for the com- 
Since all quantities calculated in lattice QCD are di- putation of the nucleon electromagnetic form factors [3. 
mensionless we need to determine a scale to convert to 



ing a is set using the nucleon mass [Tl|, [12]. The pion 
mass values, spanning a mass range from 260 MeV to 
470 MeV, are taken from Ref. (H]. At ^ 300 MeV 
and /3=3.9 we have simulations for lattices of spatial size 
L~2.1 fm and L~2.8 fm allowing to investigate finite 
size effects. Finite lattice spacing effects are studied us- 
ing three sets of results at 13=3.9, /3=4.05 and 13=4.2 for 
the lowest and largest pion mass available in this work. 
These sets of gauge ensembles allow us to estimate lat- 
tice systematics in order to produce reliable predictions 
for the nucleon form factors. 
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The authors found that cut-off effects on the nucleon 
masses are small enough to justify the application of con- 
tinuum chiral perturbation theory. Doing so the scale has 
been set through the nucleon mass at the physical point 
resulting in the following values: 

a^=3.9 = 0.089(1)(5), 
a^=4.05 = 0.070(1)(4), 
ap=A.2 = 0.056(2)(3). 

For a more detailed description, see Refs. [l^ll^. The 
mean values are systematically higher than the lattice 
spacings determined from /^r [24], but agree within one 
standard deviation. Since we are dealing with baryon 
properties we will use the values determined from the 
nucleon mass for the conversion of our lattice results to 
physical units. We note that results on the nucleon mass 
using twisted mass fermions agree with those obtained 
using other 0{a^) improved formulations for lattice spac- 
ings of about 0.1 fm and below 



III. RESULTS 

In this section we discuss the results obtained for 
the isovector electromagnetic form factors G^"{Q^) and 
(Q^), as well as the anomalous magnetic moment 
and Dirac and Pauli mean squared radii derived from 
these form factors. Before extracting values that can be 
compared to experiment we must examine the volume 
and lattice spacing dependence of these form factors. As 
already mentioned, cut-off and volume effects were stud- 
ied for the nucleon mass and taken into account in deter- 
mining the lattice spacing. We perform a similar analysis 
in the case of the form factors. 
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FIG. 



and G^j " 



2: The nucleon isovector form factors 
at rrin ~ 300 MeV for a lattice of size 24^ x 48 (filled red 
circles) and 32^ x 64 (filled blue squares). The dashed lines 
correspond to a dipole parametrization of Eq. (1221) with the 
dipole mass tue and mjv/ taken to be the p— meson mass nip 
determined on the 24'^ lattice. The dotted lines are dipole fits 
to the lattice data. The value of the magnetic form factor at 
= is fitted to the lattice data. 



A. Volume dependence 

In Fig. [2] we check for finite volume effects by com- 
paring results obtained a.t /3 — 3.9 on a lattice of spatial 
length L = 2.8 fm and L = 2.1 fm at - 300 MeV or 
for Lm,r — 3.3 and Lmj^ — 4.3, respectively. As can be 
seen, data from both volumes are compatible with each 
other for G^.j"' indicating that finite volume effects are 
are negligible. For G^"" there is an indication that the 
slope increases for the larger volume as can be seen by 
the dotted lines that are dipole fits to the form 

1 



G^ {Q )— , 2 / 2 ^2 ' 



90 



(1- 



(22) 



It is interesting to note that setting ms and tum in Eq. l22l 
to the mass of the p-meson as calculated in the lattice 
simulations yields a very good description to the lattice 
data for the larger volume. This can be seen in Fig. [21 



where the dashed line showing the dipole with the p- 
meson mass coincides with the dotted (blue) line, which 
is the fit to the data of the larger volume. The p-meson 
mass used is the one computed on the 24"^ lattice which is 
in agreement with the one computed on the larger lattice. 



B. Cut-ofF efTects 



To assess cut-off effects we compare in Figs.[3]and|3]re- 
suhs for GPr"(g2) and G57"(Q2) for three different lat- 
tice spacings at a similar pion mass. We consider results 
at our heaviest and lightest pion masses. For G5J7"(Q^)i 
results at these three lattice spacings are consistent for 
both heavy and light mass indicating that cut-off effects 
are negligible for these lattice spacings at our current 
statistical precision. There is also consistency for the re- 
sults obtained for G^~"((5^) at the smaller pion mass of 
260 MeV, as can be seen in Fig. |4l At the heavier pion 
mass of 470 MeV, a couple of data obtained on the finest 
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FIG. 3: The nucleon isovector electric (upper) and mag- 
netic (lower) form factors G^"{Q^) and GlJ"{Q'^) at ~ 
470 MeV at /3 = 3.9 (filled red circles), 4.05 (filled green trian- 
gles) and 4.2 (magenta stars) versus Q^. The open symbols 
and crosses denoted the values at — at P = 3.9, 4.05, 
and 4.2 respectively, extracted by fitting the data to a dipole 
form. 
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FIG. 4: The nucleon isovector electric (upper) and mag- 
netic (lower) form factors G|r"((5^) and G^"(Q^) at m.^ ~ 
260 MeV at /3 = 3.9 (filled red circles) and 4.2 (magenta stars) 
versus . The open symbols and crosses denoted the values 
at Q"^ = at /3 = 3.9 and 4.2 respectively, extracted by fitting 
the data to a dipole form. 



lattice have higher values, which however are well within 
the statistical fluctuations. 



C. Mass dependence of form factors 

Our lattice simulations use light quark masses that cor- 
respond to pion masses in the range of about 470 MeV to 
260 MeV. In order to obtain results at the physical point 
we need to study the dependence on the quark mass or 
equivalently on the pion mass. We show in Fig. [5] the 
dependence on the pion mass at a fixed volume and lat- 
tice spacing. We show both G^"" and G^7" computed 
at several values of the pion mass spanning pion masses 
from about 470 MeV to 300 MeV at /3 = 3.9. 

In order to extract the anomalous magnetic moment 
and mean squared radii we need to perform a fit to the 
Q^-dependence of the form factors. For both "(Q^) 
and G^J'^{Q^) we use a dipole of the form given in 
Eq. ([22]) . We fit the lattice data using all data up to 
a largest Q"^ value of ~ 1.5 GeV^. These fits, shown in 



Fig. m describe the results very well. On the same figures 
we also show Kelly's parametrization of the experimen- 
tal data [1^. Although the mass dependence is weak 
and lattice data show a weaker Q^-dependence than ex- 
perimental data, the general trend is that lattice results 
approach experiment as the light quark mass decreases 
towards its physical value. The values extracted from 
the fits for the dipole electric and magnetic masses, tue 
and TTiM 1 are therefore larger than in experiment which is 
expected given the smaller slope exhibited by the lattice 
data. As we discuss in the next section, this behavior 
is observed with other lattice discretization schemes as 
well. In Table HIl we tabulate the resulting fitting param- 
eters for all P and /i values. The parameters the, Gm(0) 
and niM have been extracted from fits to the form given 
in Eq. (j^ . 



D. Comparison of lattice results 

We have shown that volume and cut-off effects are 
small on the isovector form factors for the parameters 
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FIG. 5: The nucleon form factors G|""(Q^ 
ai P — 3.9 for = 468 MeV (crosses), m.^- 
(filled red circles) and m,r = 303 MeV (filled blue triangles) 
versus Q^. The dashed lines are the result of a dipole fit to 
the lattice data. The solid line is J. Kelly's parametrization to 
the experimental data [25j . 



used in our simulations. This justifies to some extent a 
comparison with results of other collaborations that use 
different fermions but in a similar volumes and at a sim- 
ilar lattice spacings. In Figs. [5] and [7] we show a compar- 
ison of the results of this work with those obtained using 
Np = 2 + 1 dynamical domain wall fermions (DWF) [26| , 
Np = 2 Wilson improved Clover fermions [23| and using 
a hybrid action of tadpole-improved Np = 2 -|- 1 stag- 
gered fermions and domain wall valence quarks [l^ for a 
pion mass around 300 MeV. We can see a nice agreement 
among all lattice results for G^~". As already pointed 
out, all lattice results show a weaker Q^-dependence than 
experiment. We would like to note that in Ref. (28j hy- 
brid results for — 356 MeV are obtained on a lattice 
of L = 2.5 and L = 3.5 fm. A comparison between these 
results showed volume effects for the isovector Fi but 
not for the isovector F2. This is consistent with our re- 
sults on the isovector electric and magnetic form factors 
shown in Fig. [5] The experimental data are obtained by 
interpolating the neutron form factors to the Q^-values 
of the proton form factors as described in ref. [2l|. In 
the case of G^" there are discrepancies, in particular, 
with the results using Clover fermions which are system- 



atically lower. We note that compared to G^^ " 

is more sensitive to the nucleon mass which is needed as 
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B = 


4.05 
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1.41(74) 
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1.56(19) 


2.53(49) 




0.006 


0.4035 
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4.26(30) 


1.43(20) 
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0.003 


0.2925 
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/3 = 


4.2 








0.0065 


0.4698 


1.75(12) 


3.65(24) 


2.09(31) 


2.62(24) 




0.002 


0.2622 


1.12(8) 


4.05(28) 


1.17(12) 


1.84(17) 



TABLE II: Results on the nucleon electric and magnetic mass 
extracted by fitting to a dipole form. The two last column give 
the Gm(0) and the nucleon anomalous magnetic moment in 
Bohr magnetons by fitting G^fJ"{Q^) to the dipole form given 



in Eq. ([22] 
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FIG. 6: 



Isovector electric form factor 
tion of Q^. Nf ^ 2 TMF results at 



(Q^) as a func- 
= 298 MeV are 

shown with filled squares, Nf =2 + 1 DWF [2^ (crosses), 
hybrid [2^ (open orange circles) and Clover fermions [131 (yel- 
low stars). Experimental data are shown with the filled green 
circles accompanied with Kelly's parametrization shown with 
the dashed line. 



an input. Clearly one has to study further the systemat- 
ics on the various lattice results, some of which are still 
preliminary, in order to clarify these discrepancies. 



E. Chiral extrapolation 

Given that cut-off effects are small, we use chiral per- 
turbation theory that holds in the continuum to study 
the quark mass dependence of the electromagnetic form 
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FIG. 7: Isovector magnetic form factor Gm{Q ) as a function 
of . The notation is the same as that of Fig. |6l 



factors down to the physical point. This will be justified 
in the next section where we discuss the continuum ex- 
trapolation of our results. For the chiral extrapolation, 
we use our TMF results that cover a range of pion masses 
from about 470 MeV down to about 260 MeV. The pion 
mass dependence for the isovector form factors as well 
as for the anomalous magnetic moment and radii have 
been studied within HBxPT in the so called small scale 
expansion (SSE) formulation [29| . 

The anomalous magnetic moment, which is given by 
the isovector Dirac form factor ^2(6) is extracted by fit- 
ting the assuming e.g. a dipole form dependence. 
The slope of Fi a,t — determines the transverse size 
of the hadron, < >= ~AdFi/dQ^\Q2^Q. In the non- 
relativistic limit the root mean square (r.m.s.) radius is 



related to the slope of the form factor at zero momentum 
transfer. Therefore the r.m.s. radii can be obtained from 
the values of the dipole masses by using 

(23) 

The electric and magnetic radii are given by < r^^j 
\2/mE,M and can be directly evaluated from the values 
given in Table [TTl 

Using HBxPT to one-loop, with A degrees of freedom 
and iso-vector A^-A coupling included in LO [2^, H^l the 
expression for the isovector anomalous magnetic moment 
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(24) 



and for the isovector Dirac form factor 3C 



Fr"(m.,Q2) = i 



gci - Igl - 2Bro(A) ) - 
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(25) 



To the same order the expansion of the isovector Pauli form factor is given by 



^2''^"(m„,g2) =.jf-"(m,)-g^ 
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where We perform a fit to Fi(m^,(52) i^2(™7r,Q^) with 



A /A 



k2 



R{m) ^ — + \ — -1, ^jni + Q^x{l-x). (27) 

m V m-^ 



9 



five parameters, namely the iso- vector magnetic moment 
at the chiral hmit Kv(0), the isovector and axial N to 
A coupling constants, cy and ca and the two countert- 
erms i?io(A) and E\(X). The rest of the parameters are 
fixed to their physical values, namely — 0.938 MeV, 
gA = 1.267, = 0.0924 MeV and the A-nucleon mass 
splitting A — 0.2711. The counterterms are evaluated at 
A — 0.6 GeV. For the chiral extrapolation we use data 
at the three lowest Q^-values. As can be seen in Fig. |51 
the chiral extrapolation decreases the value of Fi and 
increases the value of F2 at low Q^, bringing them into 
qualitative agreement with experiment. The values of 
the fit parameters are given in Table IIIII Although such 
a chiral extrapolation is useful and probes the general 
trend we would like to point out that the per degree 
of freedom (d.o.f) is 4.3, which means that the descrip- 
tion of the results is not really optimal. 



Fit parameter 


data at 3 /3 values 


continuum data 


Fit separately k""", r^""' 


and 


«:v(0) 
cv 

£1(0.6 GeV) 
Bio(0.6GeV) 
Bc2(0.6GeV) 


4.22(84) 
-5.86(3.52) 
-8.95(3.80) 
0.04(2) 
0.14(3) 


4.22(74) 
-5.46(2.73) 
-8.58(3.00) 
0.04(2) 
0.12(3) 


Fit to the Dirac and Pauli form factors 


«:v(0) 

CA 
Cv 

Si (0.6 GeV) 
Bio (0.6 GeV) 


5.57(56) 
1.56(3) 
-2.20(1.62) 
-6.45(2.42) 
1.21(6) 


5.20 (20) 
1.53(2) 
-3.25(53) 
-7.87(69) 
1.11(6) 
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Q2 (GeV^) 
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FIG. 8: Fi and F2 for the coarser and finer lattices at 
pion mass of about 260 MeV. The dashed line is Kelly's 
parametrization to the experimental data. The open squares 
show the values of Fi and F2 after chiral extrapolation to the 
physical point. 



TABLE III: The values of the parameters extracted from the 
fit to Fi and F2 using lattice data a.t j3 — 3.9, /3 = 4.05 and 
P — 4.2 (second column) and using the data after taking the 
continuum limit (third column). 

Using the fit parameters determined from Fi and F2 we 
can obtain the mass dependence of the isovector magnetic 
moment and radii. The expressions for the radii r^^~" 



and 



2 p—n 



are given by [3C 



1 



(47r/.)2 
A 



(l05U2)log(^ 



26-1-30 



30 



1 



logi?(TO^) 



g'A'niN 



Kv(m^) 
24r7ijv Bc2 



(28) 
(29) 



These are shown in Fig. |9] and the values obtained at the 
physical point are in agreement with experiment, which 



again indicates that the chiral extrapolation of the form 
factors could bring lattice data into agreement with ex- 
periment. 

Alternately, fitting only using three parameters, 

namely Kv(0), cy and Ei{\ = 0.6 GeV) and fixing ca — 
1.125 yields x^/d.o.f. = 0.5 and provides a nice fit to the 
results on k^^". The Dirac radius r\^~^ has only one fit 
parameter, whereas the combination r'^^~"' kP~'^ at lead- 
ing one- loop order would be predicted since the term pro- 
portional to Bc2 would be absent. However, one can allow 
for such a term, which parametrizes the short-distance 
contributions to the Pauli radius and which can be re- 
garded analogously to _Bio(A) in the Dirac radius [30| . 
We perform a fit to r^^^" with fit parameter Bio ^ind 
with fit parameter Bc2 using the expres- 



to r. 

sions of Eq. The resulting values of the parameters 

from fitting k^"", r^^~"' and r2^~^ kP~"' independently 
are given in Table IIIII If one treats ca as a fit parame- 
ter and performs a combined fit with six parameters to 
r\^~^ and r^^~^ one obtains x^/d.o.f = 1.6. 
The resulting fit for is the same as that obtained by 
fitting separately the data on k^"". The value of 
at the physical point extracted from fitting the Fi and 
F2 agrees best with the physical value. This is also true 
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FIG. 9: Chiral fits using tiie parameters determined from fit- 
ting Fl and F2. We siiow lattice data using Np ^ 2 TMF ,23] 
(a — 0.089 fm: filled red circles for L = 2.1 fm and filled blue 
squares for L = 2.8 fm; a = 0.070 fm: filled green triangles 
for L = 2.2 fm; a = 0.056 fm: purple star for L = 2.7 fm and 
open yellow square for L = 1.8 fm), iVf = 2 + 1 DWF [3l[ 
(crosses for a = 0.114 fm and L = 2.7 fm), Nf = 2 + 1 using 
DWF and staggered sea [28[ (a — 0.124 fm: open orange cir- 
cles for L = 2.5 fm and open cyan triangle for L = 3.5 fm) and 
using Nf = 2 Clover fermions [3^ are also shown with the 
cyan cross-in-square symbols. The physical point is marked 
by the asterisk. 



for the Dirac radius, whereas for the Pauli radius the de- 
viation from the physical value is larger. Omitting the 
term proportional to Bc2 in the Pauli radius increases the 
value of x^/d.o.f and the resulting fits are similar except 
for the Dirac radius above pions of about 300 MeV. 



IV. RESULTS IN THE CONTINUUM LIMIT 

In order to study the dependence on the lattice spac- 
ing quantitatively we use the simulations at three lattice 
spacings at the smallest and largest pion mass used in this 
work. We take as reference pion mass the one computed 
on the finest lattice and interpolate results at the other 
two /3- values to these two reference masses. In Fig.lTUlwe 



show the value of the Dirac and Pauli Fi and F2 at these 
reference pion masses computed in units of Tq. We note 
that we first interpolate these form factors to the same 
value of Q^. In the figure we show the form factors at 
— 0.357 GeV. We perform a fit to these data using a 
linear form Fi(a^) — i^l(O) -|- c(a/ro)^. The resulting fit 
is shown in Fig.[Tni Setting c = we obtain the constant 
line also shown in the figure. As can be seen, for both 
large and small pion masses the slope is consistent with 
zero yielding a value in the continuum limit in agreement 
with the constant fit. Therefore, we conclude that finite a 
effects are negligible and for the intermediate pion masses 
we obtain the values in the continuum by fitting our data 
a,t 13 — 3.9 and /3 = 4.05 to a constant. 















(/3 = 3.9) 


{[3 = 4.05) 


(/3 = 4.2) 


(a ^ 0) 


1.1019 


2.12(13) 


2.13(15) 


1.84(17) 


2.05(8) [1.79(26)] 


1.0 


2.08(18) 


2.36(23) 




2.18(14) 


0.95 


2.09(21) 


2.42(22) 




2.25(15) 


0.85 


2.27(27) 


2.53(26) 




2.41(19) 


0.686 


2.34(36) 


2.54(51) 




2.40(29) 


0.615 


2.49(23) 


2.57(53) 


2.62(24) 


2.55(16) [2.71(41)1 



TABLE IV: In the second, third and fourth column we give 
the interpolated values of at the value of m.^ro given in 

the first column. We used ro/a = 5.22(2), 6.61(3) and 8.31(5) 
for /3 = 3.9, 4.05 and 4.2, respectively. In the fifth column we 
give the value of after extrapolating to a = using a 

constant fit. In the parenthesis we give the corresponding 
values when using a linear fit. 

In Figs. im and . [T2]a we show the continuum extrapo- 
lation of the r.m.s Dirac and Pauli radii and the anoma- 
lous magnetic moment, respectively. The corresponding 
values of at the six reference pion masses used in 

the figures are given in Table IIVI and those of the Dirac 
and Pauli mean square radii in Table [V] 

Having results in the continuum limit we can now per- 
form the chiral fits described in the previous section. We 
show these chiral fits to the continuum results for the 
anomalous magnetic moment and Dirac and Pauli mean 
square radii in Figs. [131 HH and [151 

The behavior observed is similar to that obtained when 
using the raw lattice data. Namely, chiral fits to the Dirac 
and Pauli form factors -Fi and F2 bring agreement with 
experiment at low Q^-values, and therefore the values 
for K^'"" and " derived using the parameters of the 
chiral fit to Fi and F2 agree with the the experimen- 



tal values. The description of 



2p-n 



is also reasonable 



bringing lattice results close to the value obtained at the 
physical point, although not fully reproducing the exper- 
imental value. In the figures we also include the curves 
obtained by fitting separately the anomalous magnetic 
moment and radii. For the former the mean value ob- 
tained at the physical point is lower as compared to the 
value obtained from fitting Fi and F2, with, however, 
almost overlapping errors. For the Pauli mean squared 
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FIG. 10: Continuum results for the isovector Dirac and Pauli 
form factors Fi and F2 at = 357 GeV^. The dashed hne 
is a fit to a constant, whereas the dotted line is a fit to a line. 
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FIG. 11: Continuum results for the isovector Dirac and Pauli 
mean squared radii r\^~^ and r^^~^ . The notation is the 
same as in Fig. 1101 
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radius the fits are almost identical, whereas for the Dirac 
radius both fits do not provide a good description to our 
lattice results. The parameters extracted from the chi- 
ral fits to the continuum extrapolated results are in fact 
in agreement with those determined using lattice data 
at the different /3-values, as can be seen from the values 
given in Table Hill This is an a posteriori justification of 
using the continuum HB perturbation expressions to fit 
the lattice data at finite lattice spacing. 



FIG. 12: Continuum results for the nucleon anomalous mag- 
netic moment . The notation is the same as in Fig. 1101 
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2p-n 

3.9) 


2 p — n 


_2p-n 
^2 

4.05) 


(/3 = 4.2) 


2 p — n 

{a- 


-^0) 


1.1019 


0.236(17) 


0.300(39) 


0.229(17) 


0.319(45) 


0.183(16) 0.238(43) 


0.214(9) [0.160(39)] 


0.285(24) [0.226(70)] 


1.0 


0.226(36) 


0.379(79) 


0.258(33) 


0.326(62) 




0.243(24) 


0.347(49) 


0.95 


0.234(41) 


0.382(89) 


0.259(32) 


0.345(59) 




0.249(25) 


0.356(49) 


0.85 


0.286(30) 


0.282(81) 


0.261(35) 


0.382(68) 




0.276(23) 


0.340(52) 


0.686 


0.378(41) 


0.409(99) 


0.221(50) 


0.494(126) 




0.442(78) 


0.442(78) 


0.615 


0.307(27) 


0.446(61) 


0.210(53) 


0.534(132) 


0.266(32) 0.441(50) 


0.280(17) [0.220(54)] 


0.450(37) [0.445(91)] 



TABLE V: In the second, third and fourth raws we give the interpolated values of r^^~" and r^^'^'m ira? at the value of m.^ro 
given in the first column. We used ro/a — 5.22(2), 6.61(3) and 8.31(5) for /3 — 3.9, 4.05 and 4.2, respectively. In the fifth 
raws we give the value of r^^~" and fj^"" after extrapolating to a = using a constant fit. In the parenthesis we give the 
corresponding values when using a linear fit. 



V. CONCLUSIONS 

Computing the electromagnetic form factors of the nu- 
cleon directly from the fundamental theory of the strong 
interactions has been the goal of hadron physics since 



the discovery of QCD. Within the lattice formulation this 
goal is now being realized. Comparing lattice results with 
a number of different fermion discretization schemes we 
find an overall agreement. In this work we use dynami- 
cal simulations of two-degenerate flavors of light quarks 
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FIG. 13: Chiral extrapolation of the results on the nucleon 
anomalous magnetic moment kF^" after taking the contin- 
uum limit. The dashed line is obtained by fitting in- 
dependently, whereas the solid line is the result of fitting the 
Dirac and Pauli form factors. 
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tistical accuracy for the magnetic form factor. A small 
volume dependence is seen in the case of the electric form 
factor that indicates an increase in the slope as the vol- 
ume increases. By considering the continuum limit using 
results at three lattice spacings we also show that, cut- 
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FIG. 15; As in Fig.[T3]but for 4"-"^.^-". 



off effects are small for lattice spacings less than about 
0.1 fm. The pion mass dependence is examined using 
HB effective theory with explicit A-degrees of freedom. 
Fitting the isovector Dirac and Pauli form factors at low 
we show that the chiral extrapolated data agree with 
experiment. This is true when using the lattice results 
at the three /3-values as well as when using the lattice 
results after taking the continuum limit. 



FIG. 14: As in Fig. [13] but for 
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in the twisted mass formulation of QCD, which at max- 
imal twist is automatically O(a^) improved, thus requir- 
ing no improvement on the operator level. We use light 
quark masses yielding pion masses in the range of about 
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By comparing results at two different volumes wc find 
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